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A new search for the gravitational lens effects of a significant cosmological 
density of supermassive compact objects (SCOs) on gamma-ray bursts has 
yielded a null result. We inspected the timing data of 774 BATSE-triggered 
GRBs for evidence of millilensing: repeated peaks similar in light-curve shape 
and spectra. Our null detection leads us to conclude that, in all candidate 
universes simulated, Qsco < 0.1 is favored for 10 5 < Msco/M & < 10 9 , 
while in some universes and mass ranges the density limits £1X6 ELS much as 
10 times lower. Therefore, a cosmologically significant population of SCOs 
near globular cluster mass neither came out of the primordial universe, nor 
condensed at recombination. 
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I. INTRODUCTION 



As gamma ray bursts (GRBs) occur in the distant universe, their light is susceptible to 
gravitational lensing by intervening matter. In the gamma-ray band, the spatial resolution 
of current detectors is not great enough to allow the detection of individual GRB images, but 
the microsecond temporal resolution of many GRB detectors might allow for the detection 
of individual images temporally ffl. 

The nature of much of the dark matter in the universe remains unknown. Carr & Rees 
and Gnedin & Ostriker B suggested the theoretical possibility that condensates near the 
Jean's mass 1O 6 ' 5 M could have been the first pregalactic objects, created near the time of 
recombination in the early universe, with their mass and cosmological density depending on 
conditions. On observational grounds, globular clusters are the oldest objects known, and 
darker, more compact objects on their mass scale might well have also been created and 
escaped detection || |J. We will refer to these dark matter candidates as Supermassive 
Compact Objects (SCOs). 

Paczynski first suggested the possibility that GRBs could be seen to undergo gravita- 
tional lensing. The very short durations of some GRBs make them good sources for the lens 
detections of SCOs 0], a type of lensing known as millilensing. Previously, Nemiroff et al. 
H searched 44 GRBs for millilensing, and the null detection was used to show the universe 
was not filled to closure density of such lenses [§]. Using the same data, the (Om, ^a) = 
(0.15, 0.0) cosmology @| was marginally ruled out Marani et al. [f|. 

To act as a gravitational lens, a SCO must be compact enough to avoid absorbing one of 
the passing images. Crudely, this means that the lens must be smaller than its own Einstein 
ring. For the above canonical mass at a redshift of unity in a smooth (Qm, ^a) = (0.70, 
0.30), this radius is about 25 pc. Some massive, compact globular clusters could make good 
millilenses, but many others would not. Black holes, of course, fit all SCO and millilens 
requirements. 

In Section 2, we report the details of a new search for millilensing on 774 GRBs, a 
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factor of >15 increase in sources over the last published millilensing search ||. Again, no 
millilensing candidates were found. In Section 3, we describe the theoretical implications of 
this null detection in terms of limits on cosmological abundances of SCOs. In Section 4, we 
provide some discussions and conclusions. 



Since gravitation does not create a path-dependant time-dilation effect on passing pho- 
tons, the light-curves of different gravitational lens-induced GRB images are expected to 
be identical. Similarly, since gravitation deflects photons of different energy equally, the 
spectra of different gravitational lens-induced GRB images are expected to be identical. We 
therefore search for millilensing by comparing the spectra and the light curves of different 
peaks in the time series of GRBs. The expected time delay between GRB images created 



where z\ is the redshift of the SCO lens, / is the ratio in brightness between the brightest two 
images and is defined to be less than unity, c is the speed of light, Rs is the Schwarzschild 
radius of the SCO lens, and M is the mass of the SCO lens. Note that for canonical values 
like M = 10 7 M o , / = 0.3, and Z[ = 1.5, the expected time delay between images would be 
on the order of 20 seconds. 

The GRBs we inspected for millilensing were detected by the Burst and Transient Source 
Experiment (BATSE) onboard the Compton Gamma Ray Observatory (CGRO). This search 
was the first fully automated search. GRB data were only inspected visually if they passed 
several automated millilensing tests. Initial candidates included all GRBs detected between 
1991 April and 1999 January that had 64-ms ASCII data available on-line at the Compton 
Gamma Ray Observatory Science Support Center: 1796 GRBs in all. The last GRB in- 
vestigated was BATSE trigger number 7353. The data typically spanned about up to 524 
seconds. 



II. THE MILLILENSING SEARCH 



by a single SCO is @ 




(2.1) 
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GRBs were eliminated from the candidate millilensing sample if they had any data 
gaps. GRBs with data gaps would likely not show any type of millilensing bias but would 
pose special challenges for a program written in a general, automated fashion. Surviving 
GRBs were then eliminated from the candidate millilensing sample if they had backgrounds 
that sloped by more than 4.5 a between the first 64-ms time bin and the last 64-ms time 
bin. Again, these GRBs should show little bias toward millilensing but posed unwieldy 
programming challenges. After these cuts 774 GRB millilensing candidates survived. 

The time series of all surviving GRBs were then analyzed by an automated millilensing 
search program: millisearch.for. The program's first task was to isolate the first discrete 
episode of GRB emission, which we will call the "burst." To start, the highest 64-ms time 
bin of each surviving GRB was first found. This high bin was then added to the very 
next trailing bin to see if the signal to noise (S/N) above background was increased. If 
not, the burst was considered to have ended. If so, the trailing bin was considered part of 
the burst, and added to the burst. The next trailing bin was then considered and treated 
identically. After the last trailing bin was added to the burst, leading bins were then tested 
and added, again until the S/N over background was found to decrease. At this point, the 
first contiguous episode of GRB emission was taken to be defined. 

Next, the program searched to see if a second identical (to within an intensity scale 
factor), gravitational- lens induced image of the burst could possibly have followed in the 
available BATSE time series. A candidate emission episode is referred to here as an "echo." 
To find an echo, the highest 64-ms time-bin following at least one full burst duration after the 
conclusion of the burst was first isolated. A series of 64-ms time bins around the echo peak 
was then isolated equal to the duration of the burst. GRBs where the echo was below 5.5 a 
over the background were recorded as having no significant echo. For those GRBs without a 
candidate echo, maximum echo strength f max , minimum time delay between burst and echo 
t m in, and maximum time delay between burst and echo t max were recorded for the analysis 
described in the next section. 

Those GRBs with a candidate echo had f max recorded as the counts ratio between the 
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echo and the burst. As before t m i n and t max were also recorded for these GRBs. All the 
counts in all the time bins of the burst and echo were then co-added in all four major BATSE 
energy channels. The relative counts in each energy channel were compared between burst 
and echo by a two-distribution \ 2 test- Eleven GRBs passed this test, and were then 
compared for light-curve similarity. One such candidate is shown in Figure 1. No GRBs 
were found to have both matching millilensing spectra and light-curves. The final output of 
the millilensing search program was a list of BATSE trigger numbers, and f max , t m i n , and 
t max for each trigger. 

Faux GRBs with a fictitious millilensing signal were used to successfully test the millilens- 
ing detection program. This program, millisearch.for, is available through the Astrophysics 



Source Code Library at |http:/ /ascl.net 



III. COSMOGONIC IMPLICATIONS 

As millilensing would be expected in a variety of cosmogonic scenarios, non-detection 
rules out some of these scenarios, and sets limits on the potential abundance of SCOs in 



other scenarios. To compute these limits, we used the detection volume formalism [|TT|] ||. 
To be detectable, lenses must fall inside a volume where they would create two bright images 
with dynamic range between them f m i n . Additionally, lenses must fall inside another volume 
where the two bright images would be separated in time by less than a maximum amount 
tmax, and fall outside a third volume where the two bright images would be separated in 
time by less than a minimum amount t m i n . 

At each lens redshift zi, a detection volume has a radius b. For the detection criterion 



of magnification, the detection volume radius is given in JTTJ as 6 a = J ARsDq L D^ s ^ / Dps 
where Rs is the Schwarzschild radius of the compact lens, D A refers to angular diameter 
distance, and subscripts O, L, SzS refer to the observer, lens and source respectively. Magni- 
fication A is embedded in $ = A/\/ A 2 — 1 — 1. Dynamic range / between the two brightest 
images of a compact lens is by definition less than unity and related to total magnification 
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by A = (1 + /)/(l — /). Similarly, the time delay At between the two brightest images of a 
compact lens is related to dynamic range by Eqn. 2.1. 

A lens is detectable if it falls inside both bf and b tma x but outside btmin- The total number 
of detectable millilenses is therefore given by 

Ni ens = niir(b 2 uter - b\ nner )dD p OL (3.1) 

J 

where Dq L is the proper distance between the observer and the lens, ni is the comoving lens 
density, b oute r is the radius of the outermost detection volume where b oute r = min (bf, btmax) 
and bi nner btmin- 

In a few cases, GRB source redshifts z s have been actually determined (see, for example, 



[ 12 1 ) . To find the likely effects of millilensing, however, we needed redshifts for the whole 
BATSE sample. We therefore estimated these redshifts from Monte Carlo simulations. These 
simulations are described in detail elsewhere ||14|| , and have the following characteristics. 



GRBs are thrown randomly in a universe with a specified geometry. These GRBs are 
assigned a generic spectrum and a peak flux chosen from a published luminosity function 
[ 375 1 . These GRBs are then redshifted and re-sampled to see if they would re-trigger BATSE. 



Those that do are fitted to the BATSE 3B detection rate, the BATSE Log N - Log P, time- 
dilation [17 , and (V/V max ). The later statistic supplemented the normalization of the 



Log N - Log P fit, which was most sensitive to brightness distribution shape. Adjustable free 
parameters included an integrated flux scale factor L Q and a scale factor on the comoving 
number density n s . A file recording acceptable fits for each universe geometry was recorded. 
Each of the 774 BATSE millilensing candidates was then matched-up to the Monte-Carlo 
GRB with the closest peak flux, and assigned the redshift of this Monte-Carlo GRB. 

Results are shown for two of these cosmologies in Figures 2 and 3. In the figures, the 
number of expected lens detections is plotted as a function of lens mass and density. The 
contours represent 1, 5, and 10 expected lens detections, and hence roughly correspond to 
exclusions at the 1, 2, and 3 a levels. 

The expected number of lens detections in a (Qm, Q\) = (0.3, 0.7) universe |1| [19| is 
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shown in Figure 2, on the left. VLm is the density of matter in the universe, relative to the 
critical density, while Vt^ is the density in a cosmological constant. 

Here one finds a wide range of masses is excluded for several interesting values of ttco- in 
fact, Vlco dips below even 0.03 for a narrow range of masses near M/M & ~ 10 7 . Therefore, 
even a significant amount of baryonic dark matter is excluded from forming compact lenses 
in this mass range. 

To test Gnedin & Ostriker's || universe directly, Figure 3, on the right, shows similar 
results for a (0.15, 0) universe. It appears difficult for any significant fraction of the mass 
density to hide in supermassive compact objects in the mass range 10 6 < M/M Q < 10 8 . 



IV. DISCUSSION AND CONCLUSIONS 

Lensing of GRBs has the ability to time-resolve a mass range of dark matter not usually 
probed: that of compact objects near globular cluster mass. That these lenses were not 
detected indicates, as we have shown above, that the majority of the universe's mass density 
does not reside in SCOs. 

A potential inaccuracy in the above analysis lies in the estimated GRB redshifts. The 
average GRB redshift we have estimated is between 1.5 and 2, like most so far determined in 
the optical. However, some popular speculation holds that average BATSE redshifts would 
be higher, and, if true, would make the limits estimated here even more conservative. 

The above lens rates did not incorporate the potential effects of magnification bias. 
Such bias [^fj] would be expected were the GRB luminosity function to rise sharply near the 
effective redshift limit of the BATSE sample. One reason for not including magnification 
bias is the potential uncertainty in the slope of the GRB luminosity function. As including 
magnification bias might significantly increase the expected amount of lensing, it makes the 
estimated limits derived here even more conservative. 

Another potential inaccuracy would occur were different GRBs images to have different 
light-curves or spectra. This might occur were GRBs to undergo beaming [T2"T|-|2"3"|] on the 
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order of milli-arcseconds or less, something that has not been ruled out by observation. 

As more GRBs become available for millilensing inspection, the ability to find a smaller 
and smaller abundance of SCOs becomes possible. Potentially, as much as a factor of 10 
more GRBs might become available, allowing an abundance as low as flsco ~ 0.003 to be 
discernable. A survey sensitive to more GRBs and greater time resolution might detect the 
known star-field in GRB lensing, as well as dark matter on scales near one solar mass Q. 
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Note added in proof: After most of this work was completed, we became aware of work 
by Wilkinson and collaborators who were making a similar investigation into cosmological 
millilensing but using observations with VLBA radio telescopes. That group reports their 
work and results in the next paper. 
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FIGURES 



11 



Counts Measured 
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FIG. 1. BATSE GRB trigger 5457 light-curve showing two emission episodes that comprised 
one of the better candidates for millilensing. The counts in four energy channels are shown as a 
function of time. The time resolution of the data and plot is 64-ms, except before trigger, where 
it is 1024-ms. 




FIG. 2. The number of expected millilenses in a (0.3, 0.7) universe as a function of millilens 
mass, density, and universe geometry. 
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FIG. 3. The number of expected millilenses in a (0.15, 0) universe as a function of millilens 
mass, density, and universe geometry. 
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